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ABSTRACT 


4/ 

A mathematical  description  of  ambient  sea  noise  is  presented  (hat 
takes  into  account  the  two  dominant  sources  — surface  agitation  and 
distant  shipping.  The  contribution  from  surface  agitation  is  shown  to  be 
linearly  dependent  upon  the  logarithm  of  wind  speed  for  depths  between 
400  and  2500  fm  at  several  sites  near  Bermuda.  When  this  is  taken  into 
account,  the  ambient- noise  data  that  include  both  sources  can  be 
analyzed  to  determine  the  individual  levels  of  the  sources.  Also,  the 
standard  deviations  of  the  two  sources  can  be  determined  from  the 
standard  deviations  of  the  measured  levels.  r J T5t§TRlfeUTl5N  &TAT1 
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INTRODUCTION  AND  SUMMARY 

Although  many  sources  of  ambient  sea  noise  have  been  identified,  in  most 

areas  two  sources  dominate  between  10  and  10,000  Hz.  One  source  is  non-wind- 

dependent  and  believed  to  be  primarily  distant  shipping.  The  other  source  is  wind- 

dependent  and  believed  to  be  surface  agitation  near  the  measuring  device.  This  is 

borne  out  in  a survey  of  many  ocean  areas  by  Wenz  and  in  measurements  made  near 
2 3 

Bermuda  by  Perrone.  ' 

Piggott,  using  data  taken  in  very  shallow  water,  found  that  the  noise  spectrum  level 
attributed  to  the  wind-dependent  source  was  proportional  to  the  logarithm  of  wind  speed.4 
His  study  was  conducted  on  the  Scotian  Shelf  in  20  and  28  fm  cf  water.  However, 
because  of  Hie  greater  complexity  of  deep-water  data,  his  findings  have  not  been 
previously  substantiated  for  other  locations  and  other  depths.  We  have  found  that 
independently  obtained  sets  of  data  taken  at  several  sites  near  Bermuda  in  depths  up 
to  2500  fm  also  demonstrate  the  logarithmic  relationship. 


The  conclusions  that  only  two  sources  are  important  and  that  the  contribution 
from  the  wind-dependent  source  is  linearly  related  to  the  logarilhm  of  wind  speed 
simplify  ambient-noise  prediction  and  analysis.  The  former  conclusion  has  been 
taken  into  account  for  some  time,  but  it  still  liar,  been  difficult  to  separate  the 
relative  contributions  of  the  two  sources  in  measured  levels.  And  it  has  been  equally 
difficult  to  estimate  the  standard  deviations  of  the  two  contributions  from  the  standard 
deviations  of  the  measured  levels.  However,  the  linear  dependence  that  Figgott^ 
discovered,  and  which  is  substantiated  here,  makes  it  possible  to  express  mathemati- 
cally the  noise  levels  end  standard  deviations  as  functions  of  wind  speed  for  a given 
frequency.  The  spectrum  levels  and  the  standard  deviations  of  Hie  contributions  from 
both  sources  can  be  obtained  by  finding  the  mathematical  equation  of  the  appropriate 
form  that  best  fits  the  data. 
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Section  I provides  ihe  necessary  mathematical  expressions.  Section  II  demon- 
strates that  the  linear  relationship  does  exist  in  several  sets  of  deep- water  data. 

And  Section  III  discusses  applications  of  the  mathematical  relationships  in  determining 
the  spectrum  levels  and  the  standard  deviations  of  tlu  contributions  from  the  two  sources. 


I.  THEORY 

Piggott  studied  the  wind  dependence  of  ambient  sea  noise  at  two  depths  in 
very  shallow  water  (20  and  20  fin)  where  the  contribution  from  the  wind-dependent 
source  was  dominant  at  nearly  all  wind  speeds  and  all  frequencies.^  He  found  that 
at  wind  speeds  where  the  contribution  from  the  wind-dependent  source  was  dominant, 

the  spectrum  level  (in  decibels)  was  linearly  dependent  on  the  logarithm  of  wind 
speed  (Eq.  2,  below).  His  data  also  suggest  that  the  measured  levels  at  all  wind 
speeds  and  at  all  frequencies  could  be  described  as  the  sum  of  the  contributions 
from  two  sources  — the  wind-dependent  source  and  the  non-wind-dependent  source. 

Mathematically,  measured  levels  can  be  represented  as  the  sum  of  two 
portions  (Ihe  notation  is  chosen  to  be  consistent  with  Piggott's  notation.^): 

L,(f)  = A(f)  (1) 

L2(f,  v)  = R(f)  + 20  n(f)  log  v,  (2) 


where 


L,(f) 


spectrum  level  in  decibels  of  the  non-wind-dependent  portion 
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L2(f,  v) 


spectrum  level  in  decibels  of  the  wind-dependent  portion 
(as  published  by  Piggott) 


A(f)  spectrum  level  in  decibels  of  the  non-wind-dependent  portion 


B(f) 

20n(f) 


spectrum  level  in  decibels  of  Hie  wind -dependent  portion 
at  1 knot  wind  speed 

slope  of  the  wind-dependent  portion  in  decibels  per  decade 
wind  speed  in  knots 


The  level  resulting  from  random  phase  summation  of  Lj  and  I_2  is 


L(f/  v)  = 10  log 


(ioL,/,0+  io'27'0)  . 


(3) 


II.  ANALYSIS 


In  Piggott's  study,  the  wind-dependent  portion  was  dominant  at  wind  speeds 
above  25  knots  for  frequencies  from  8 to  140  Hz,*  at  higher  frequencies  (up  to  3000  Hz), 
the  wind-dependent  portion  was  dominant  at  all  wind  speeds.  He  was  able  to  fit 
straight  lines  to  the  wind -dependent  portions  by  eye  and  thus  confirm  that  the  linear 
relationship  existed.  One  would  expect  the  linear  relationship  he  found  to  hold  at 
other  locations,  but  it  would  not  be  so  pronounced  if  the  non-wind- dependent 
portion  was  relatively  higher. 


Data  have  been  taken  by  New  London  Laboratory  scientists  with  several  arrays 

2 3 

near  Bermuda  that  show  ambient-noise  levels  as  a function  of  wind  speed.  Perrone  ' 
provided  the  most  extensive  data,  and  data  also  were  available  from  two  other  studies. 
The  hydrophones  were  at  depths  between  400  and  2500  fm. 


j 
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In  oil  ca.es,  when  these  data  were  plotted  as  spectrum  levels  in  decibels  vs 
logarithm  of  wind  speed,  it  appeared  that  the  wind-dependent  portion  might  exhibit 
a linear  dependence  at  higher  wind  speeds  for  some  frequencies.  If  was  not  possible, 
however,  to  choose  the  appropriate  straight  line  by  eye  as  Piggotf4  had  done.  At 
many  frequencies,  the  wind  dependence  consisted  only  of  a small  rise  in  level  at  the 
highest  wind  speeds. 

In  order  to  find  the  parameters  B(f)  and  n(f)  of  the  linear  relations!) ip  for  these 
data,  we  used  the  mathematical  relationship  of  Eq.  3.  The  appropriate  B(f)  and  n(f) 
can  be  determined  by  finding  a best-fit  curve  of  that  form.  The  advantage  of  this 
method  is  that  if  uses  all  data  points  even  though  many  of  them  represent  the  pre- 
dominantly non-wind -dependent  portion.  The  non-wind-dependent  portion  L^f)  also 
is  calculated  in  the  analysis.  We  defined  the  best-fit  curve  as  that  curve  for  which 
the  sum  of  the  squares  of  the  vertical  distances  between  the  data  points  and  the 
curve  was  minimized. 

At  all  frequencies  where  there  was  evidence  of  a wind-dependent  portion  • 

(generally  at  10  to  15  Hz,  25  to  40  Hz,  and  above  100  Hz),  such  a curve  could 'be 

2 

fitted  to  all  the  sets  of  data.  Perrone's  data  are  used  here  to  demonstrate  the  results. 

The  data  and  the  best-fit  curves  are  shown  in  Fig.  1.  Each  curve  is  defined  by  the 
three  parameters  A(f),  3(f),  and  n(f),  and  a different  set  applies  to  each  frequency 
(see  Table  I). 

Surprisingly  little  wind  dependence  was  needed  for  a best-fit  curve  to  be 
calculated  showing  the  character  of  the  wind  dependence.  Only  when  there  was 
no  increase  in  the  spectrum  levels,  even  at  the  highest  wind  speeds,  was  it  impos- 
sible to  get  a fit.  However,  when  there  was  minimal  wind  dependence,  B(f)  and 
n(f)  were  not  reliable.  From  Perrone's  data,  2 reliable  B(f)'s  and  n(f)'s  were 
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obtained  for  11  Hz,  28  Hz,  and  for  all  frequencies  above  112  Hz.  At  14,  28, 
and  89  Hz,  unreliable  B(f)'s  and  n(f)'s  v/ere  obtained;  at  the  other  frequencies,  the 
non-wind -dependent  portion  was  dominant  even  at  a 50-knot  wind  speed.  Similar 
results  were  obtained  with  the  other  sets  of  data. 

The  data  generally  fall  at  small  random  distances  from  the  best-fit  curves; 
however,  two  small,  but  systematic,  variations  can  be  noted.  One  could  be 
described  as  a "waviness"  of  the  data.  At  the  higher  wind  speeds,  the  data 
for  many  frequencies  meander  about  the  best-fit  curve  in  the  same  manner.  This  is 
evident  at  1414  and  2816  Hz  in  Fig.  1.  The  variation  is  only  1 or  2 dB  and  may  be 
the  result  of  data-sorting  and  editing  techniques.  However,  a physical  explanation 
can  not  be  ruled  out. 

Also,  there  is  some  indication  that  the  data  may  flatten  out  at  the  higher  wind 

speeds.  This  characteristic  can  be  seen  in  the  data  for  707  Hz  in  Fig.  1.  Piggott's 

4 

results  do  not  provide  any  information  for  v/ind  speeds  above  35  knots.  The  apparent 
flattening  appears  above  that  wind  speed.  A possible  explanation  cf  the  flattening  is 
that  a maximum  level  of  the  wind-dependent  portion  is  reached,  but  the  effect  is  not 
pronounced  enough  to  warrant  any  conclusions  at  this  time. 

FJgure  2 shows  n(f)  as  a function  of  frequency  for  all  the  data  analyzed.  The 
sets  of  data  are  labeled  by  depth:  400,  1 100,  and  2500  fm  are  from  an  article  by 
Perrone,  700  to  1 100  fm  and  2100  fm  are  from  unpublished  reports,  and  2400  fm  is 

2 

from  the  set  of  data  by  Perrone  that  is  used  in  this  article  to  demonstrate  the  results. 
A function  represented  by  two  line  segments  seems  to  be  suggested  by  the  data  for-the 
four  greatest  depths.  The  line  segments  correspond  to  two  of  the  three  derived  by 

4 jr 

Piggott,  but  they  are  lower  (see  Fig.  3).  Also  shown  are  n(f)'s  found  by  Payne  in 
20  fm  of  water. 
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Froin^the  six  sots  of  data  analyzed  in  this  study  and  the  two  sets  of  data  previously 
published,  ' one  can  see  that  n consistently  falls  into  an  orderly  pattern  as  a 
function  of  frequency.  In  most  cases.  Hie  pattern  is  similar  to  that  originally  determined 
by  Piggott  in  20  and  20  fm  of  wator.  However,  the  pattern  is  found  to  be  lower  in 
Puyne's  data  for  20  fm  of  water  and  in  four  sets  of  New  London  Laboratory  data  (1 100 
through  2500  fm).  Also,  two  sets  of  New  London  Laboratory  data  (-100  and  700  Hirough 
1 100  fm)  do  not  conform  well  to  the  shape  of  the  pattern.  There  is  no  obvious  depth 
effect,  nor  is  it  otherwise  clear  why  the  differences  exist.  The  relatively  large  values 

of  n(f)  around  89  1 Hz,  which  appear  in  several  sets  of  data,  are  also  curious  and  not 
yet  explained. 


III.  applications 

A.  S pectrum  Levels  of  the  Two  Sources 

When  measured  ambient-noise  levels  are  represented  by  Lq.  3,  Lj(f)  and 
L2(f,  v)  ate  estimates  of  the  two  component  portions.  In  Section  II,  we  showed  that 
Perrone's  data  could  be  so  represented,  and  we  calculated  the  parameters  A(f),  B(f), 
and  n(f),  which  arc  necessary  for  specifying  L,(f)  and  L2(f,  v),  and,  thus,  L(f,  v) 
for  each  frequency. 

Figure  4 shows  the  calculated  levels  of  the  non-wind-dependent  portion,  L^f), 
and  Hie  wind-dependent  portion,  L^f,  v).  The  latter  portion,  given  by  Eq.  2 using 
B(f)  and  n(f)  from  Table  1,  is  plotted  for  wind  speeds  of  1,  5,  10,  20,  .30,  40,  and 
50  knots.  Although  B(f)'s  and  n(f)'s  were  found  for  14,  28,  and  89  Hz,  they  were 
judged  unreliable  and  are  not  used  here.  Calculations  wore  made  for  1 1 and  35  Hz 
and  112  through  1414  Hz,  and  L2(f,  v)  was  estimated  by  extrapolation  at  the  other 
frequencies.  The  non- wind- dependent  portion  is  given  by  Eq.  1 using  A(f) 
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from  Table  I.  One  can  readily  see  that  when  the  wind-  and  non-wind-dependent 
portions  are  combined  they  give  Hie  measured  levels. 

The  process  of  "least-square"  fitting  a curve  of  the  shape  given  by  Eq.  3 is 
mathematically  complex.  Therefore,  simpl ified  techniques  were  considered.  In 
extreme  cases  when  all  data  at  a given  frequency  are  wind -dependent,  as  in 
shallow-water  data  (Piggott\  or  non-wind-dependent,  as  it  is  at  several  frequencies 
in  deep  water,  a suitable  estimate  of  the  dominant  source  may  be  obtained  "by  eye" 
or  by  simple  straight-line  fitting  techniques.  In  other  cases,  where  contributions  of 
both  sources  are  significant,  the  noise  level  for  wind  speeds  below  about  5 knots  will 
generally  be  predominantly  non-wind-dependent,  and  may  be  taken  as  a rough 
estimate  of  Lj  in  Eq.  3.  This  suggests  that  a rough  estimate  of  can  be  obtained 
by  (1)  subtracting  from  each  data  point  and  then  (2)  performing  a simple  straight- 
line  fit.  Unfortunately,  such  an  analysis  is  unsatisfactory.  Considerably  less  data 
are  used  in  the  estimate  of  L^,  which  means  greater  probable  error,  and  errors  in  Lj 
are  amplified  at  low  wind  speeds  when  it  is  subtracted  from  the  data.  The  resultant 
error  in  proves  unacceptable. 

B.  Standard  Deviations  of  the  Two  Sources 

When  Eq.  3 applies,  the  variation  of  measured  levels  about  their  mean  would 
result  from  the  combined  variations  of  the  two  portions.  A mathematical  model  for 
the  combination  is  given  below  in  terms  of  the  standard  deviations  of  two  Gaussian 
distributions. 

The  standard  deviation  of  the  measured  levels  is  given  by 


where 
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L’  + 

,ow  ,ow 


Lj  , Lj,  and  L arc  given  by  Eqs.  1,  2,  and  3.  This  is  only  an  approximate 
relationship  but  adequate  for  the  data  currently  available. 

Figures  5 and  6 are  the  results  of  applying  the  above  model  to  the  standard 

2 # ■ 

deviations  associated  with  Perronc's  data  for  all  frequencies  and  for  wind  speeds  of  5, 
20,  30,  40,  and  50  knots.  Figure  5a  shows  the  measured  standard  deviations. 

Since  it  is  not  possible  to  solve  explicitly  for  sets  of  and  a ^ / we  first 
estimated  then  calculated  (Tj  . We  did  this  for  frequencies  at  and  above 

224  Hz  because  the  model  can  not  agree  with  the  data  at  frequencies  where  the 
measured  levels  were  totally  non-wind  dependent.  The  data  show  a wind  dependence 
in  tire  standard  deviation  of  those  frequencies,  but  the  model  shows  ° - <r^(f)  . 

Figure  5b  illustrates  the  estimate  of  (for  224  llz  and  above)  that  gave  the  most 
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consistent  results  for  a.^  - The  tr,'s  below  224  Hz  were  chosen  os  cquol  to  the  measured 
values  of  o for  5-hnot  wind  speed.  The  ° 2%  ca*cijla,ecl  are  shown  in  Fig.  5c. 
Perronc  concluded  that  was  probably  independent  of  frequency  we  found 
no  consistent  dependence.  Thus,  the  dashed  line  in  Fig.  5c  is  taken  as  the  best 
estimate  of  <J . 


2 * 


The  validity  of  Hie  analysis  is  checked  by  combining  the  a ^'s  from  Fig.  5b 


and  the  a^i  from  Fig.  5c.  Figure  6 shows  the  <r's  thus  calculated.  They  are 


a good  approximation  of  the  measured  values.  The  only  major  difference  is  the 
non-wind-dcpendcncc  of  <r  at  frequencies  where  spectrum  levels  were  also  found 
to  be  non-wind-dependent. 


This  method  of  determining  a ^ and  a ^ 's  not  exact,  and  one  is  tempted  to 


shift  a | and  o ^ slightly  to  get  a better  fit.  We  did  not  find  this  to  be  useful 


because  of  the  limited  data  available  to  us.  Researchers  with  data  that  have  not 
been  smoothed  and  that  cover  a wider  frequency  range  may  have  more  success  in 
determining  reliable  estimates  of  and  tr^  • The  primary  purpose  of  this 
discussion  is  to  demonstrate  the  usefulness  of  the  technique.  The  results  should  be 
taken  only  as  a tentative  estimate  of  what  and  a ^ may  be. 
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Table  I 


PARAMETERS  OF  THE  BEST-FIT  CURVES 


Frequency 

(Hz) 

> 

3 

% 

B(0 

n(0 

11 

-11.31 

-64.35 

1.73 

14 

-15.05 

-81.23 

2.07 

17 

- 7.01 

22 

- 8.42 

28 

-16.48 

-88. 53 

2.09 

35 

-18.59 

-69.87 

1.38 

44 

-17.  73 

56 

-14.03 

70 

-17.49 

89 

-27.69 

-125.86 

2.77 

112 

-28. 03 

-73.  17 

1.30 

141 

-30. 58 

-74.  80 

1.39 

177 

-34. 20 

-55. 54 

.81 

224 

-37.  14 

-59.01 

. 93 

281 

-38.87 

-57.64 

.87 

354 

-38.61 

-58.  1 1 

.93 

446 

-41.36 

-56.93 

.90 

562 

-44.27 

-58.  72 

.92 

707 

-49.49 

-65.66 

1.09 

891 

-50.  77 

-69.  82 

1.  15 

1122 

-52.51 

-69.  25 

1.03 

1414 

-55. 82 

-69.43 

.97 

1778 

-54.54 

-71.  15 

.96 

2241 

-52.  99 

-71.41 

.93 

2816 

-54.  79 

-73.95 

.96 

v i 
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Laboratory  of  the  Naval  Underwater  Systems  Center. 
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Fig.  1.  Ambient  noise  soectrum  levels  vs  the  logarithm  of  wind  speed  showing  the 
best-fit  curves.  The  non-wind-dependent  and  wind-dependent  portions 
of  the  best-fit  curves  are  shown  by  the  dashed  lines.  I 
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Fig.  3.  Representative  slopes  of  the  wind-dependent  portion  for  the  New  London 
Laboratory  data  compared  with  the  slopes  found  by  Piggott  and  Payne. 
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Fig.  4.  Calculated  spectrum  levels  of  the  non-wind-dependent  and  wind 
dependent  portions  as  derived  from  measured  levels. 


